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Abstract. We consider time-dependent Schrodinger equations in one dimen- 
sion with double well potential and an external nonlinear perturbation. If the 
initial state belongs to the eigenspace spanned by the eigenvectors associated 
to the two lowest eigenvalues then, in the semiclassical limit, we show that the 
reduction of the time-dependent equation to a 2-mode equation gives the dom- 
inant term of the solution with a precise estimate of the error. By means of 
this stability result we are able to prove the destruction of the beating motion 
for large enough nonlinearity. 



1. Introduction 

Recently, the theoretical analysis of the nonlinear time-dependent Schrodinger 
equation 

ifiV' = ifoV' + e|V'lV, ee R, ^ = ^, (1) 

where 

has attracted an increasing interest (see for a review and for a rigorous 
derivation of the Gross-Pitaevskii energy functional). When y is a double well 
potential, one of the main goal is to understand how the nonlinear perturbation with 
strength e affects the unperturbed beating motion (see, e.g., the review paper 
and the paper Q where equation (|l]) is proposed as a model for chiral molecules). 
To this end, it is crucial to study the solution for times of the order of the beating 
period; in other words, for practical purposes the unit of time is given by the beating 
period T ~ iih/Lj where h is the Planck's constant and uj is one half of the energy 
splitting between the two lowest energies. 

Here, I consider equation (|l|) in the semiclassical limit where, by assuming that 
d — 1 and under some generic assumption on the double well potential, we give 
the asymptotic behavior of the solution ■0 with a precise estimate of the error. In 
particular, the main result (Theorem 3) consists in proving that the solution of 
the Gross-Pitaevskii equation is approximated, with a rigorous control of the error. 
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by means of the solution of a two-dimensional dynamical system exactly solvable. 
As a result it follows (Theorem 4) that the beating motion between the two wells 
of a state initially prepared on the two lowest eigenstates gradually disappears for 
increasing nonlinearity. 

A similar investigation has been recently performed in where the nonlinear 
perturbation is given by e{ip, gip)gip and g{x) is a given odd function, and in JlSf , 
where, in dimension d = 1 and d = 3, we consider the limit of large barrier between 
the two wells. In particular, in [ p^ I had to assume that the discrete spectrum of 
the Schrodinger operator Hq consists of only two non-degenerate eigenvalues and 
that the restriction to the continuous eigenspace of the unitary evolution operator 
satisfies to a priori estimate uniformly with respect to the parameters of the model. 

Finally, we mention other recent results concerning the study of the existence of 
stationary solutions for Gross-Pitaevskii equations with double well potentials 
(2I and, in the case of single-well type potentials, the existence of solutions asymp- 
totically given by solitary wave functions in the case that the discrete spectrum of 
the linear Schrodinger operator has only one non-degenerate eigenvalue [p^ , . 

Our paper is organized as follows. 

In Section 2 we introduce the main notations and we state the assumptions 
on the potential. Moreover, we collect some semiclassical results concerning the 
spectrum of the linear Schrodinger operator. 

In Section 3 we prove the global existence of the solution of the Gross-Pitaevskii 
equation, the existence of conservation laws and a priori estimate (Theorem 2) . The 
global existence of the solution is proved for both repulsive and attractive nonlinear 
perturbation, where, in the second case, we have to assume that the strength of the 
nonlinear perturbation is small enough. 

In Section 4 we introduce the two-level approximation which, roughly speaking, 
consists in projecting the Gross-Pitaevskii equation onto the two-dimensional space 
spanned by the eigenvectors of the linear Schrodinger operator associated to the 
two lowest eigenvalues. For practical purposes, it is more convenient to choose, as a 
basis of such a two-dimensional space, the two "single- well" states. The dynamical 
system we obtain is exactly solvable. 

In Section 5 we give our main result (Theorem 3) proving the stability of the 
two-level approximation. Here, we make use of the comparison criterion between 
ordinary differential equations and of a priori estimate of the solution of the Gross- 
Pitaevskii equation. We underline that, in order to obtain such an estimate, the 
assumption d = 1 on the dimension plays a crucial role. 

In Section 6 we give the full rigorous justification of the results by Vardi [|6| 
proving the existence of a critical value for the nonlinearity parameter giving the 
destruction of the beating motion (Theorem 4). 

2. Assumptions and preliminary results 
Here, we consider the Cauchy problem 

ih^ = i/eV', H, = Ha + W (2) 

where ip denotes the derivative of V' with respect to the time t, Hq is the linear 
Schrodinger operator formally given by (here, x denotes the spatial variable in 
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dimension 1) 



y is a symmetric double well potential and 

is the nonlinear perturbation with strength e. 

In the following, for the sake of definiteness, we denote by C any positive con- 
stant independent of e, h and we assume h small enough, that is ?i G (0, h*] for 
some h* ^ and we denote 

M\p = Ml. = WixWdx^ and = H^l^. 

Moreover, given y = {yi, . . . , y^) G R™ for some m > 1, we denote 

\y\ = max \yj\. (4) 

l<j<m 

2.1. Assumptions on the potential. Here, we assume that the potential V is 
a regular symmetric function which admits two non-degenerate minima and it is 
bounded from below. More precisely: 

Hypothesis 1. The potential V{x) is a real valued function such that: 

i. V{-x) = V{x), Vx G R; 
ii. V G C2(R); 

Hi. V{x) admits two non- degenerate minima at x ^ ±a, for some a > 0, such 
that 

V{x) > V,mn = V{±a), Vx eR, x^ ±a; (5) 
in particular, for the sake of definiteness, we assume that 

^ = and -^5-^ > 0; 

dx dx'^ 

iv. finally we assume that 

liminf|a.|^oo't^(2;) = Kx) > Vmtn- 

It follows that the operator formally defined in (||) admits a self-adjoint realiza- 
tion (still denoted by Hq) on i^(R) (see, for instance, Theorem III. 1.1 in |^). Let 
cr{Ho) — CTdUCTess be the spectrum of the self-adjoint operator Hq where ad denotes 
the discrete spectrum and a ess denotes the essential spectrum. From Hyp. 1-iv it 
follows that ad C (Knm, Voo), (^ess = if = +00 (see Theorem XIII. 67 in 
and that a^ss ^ Woo^ -foo) if Voo < oo (see Theorem HI. 3.1 in |^). Furthermore, 
the following two Lemmas hold: 

Lemma 1. Let ad be the discrete spectrum of Hq- Then, for any h G (0, ft*] it 
follows that: 

i. ad is not empty and, in particular, it contains two eigenvalues at least; 
ii. let Ai,2 be the lowest two eigenvalues of Hq, they are non-degenerate, in 
particular Ai < A2 and there exists C > 0, independent of h, such that 

inf [A - A2I > Ch. 
xeT{Ho)-{Xi.2} 
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Proof. The proof is an immediate consequence of the above assumptions and stan- 
dard WKB arguments. □ 

Lemma 2. Let ^\,2 be the normalized eigenvectors associated to Ai_2, then: 

i. ipj, j = 1,2, can be chosen to be real-valued functions such that ipj{'-x) = 

a. ipj e i?^(R); 

Hi. ipj g LP(R) for any p € [1, +c»]; 

iv. there exists a positive constant C such that 

\\fj\\p<Ch-~'^, Vpe[2,+oo], vne(0,r]. (6) 

Proof. Property i. immediately follows from assumption Hyp. 1-i. Property ii. 
follows from Lemma IIL3.1 in Property Hi. follows from Theorem in.3.2 in Q. 
Finally, property iv. follows for p = +00 by means of standard WKB arguments. 
From this fact, from the normalization of the eigenvectors and from the Holder 
inequality then property iv. follows for any p G [2, +00]: 



\'Pj\ 



2 p-2| 



□ 



2.2. Splitting and single-well states. It is well know that the splitting be- 
tween the two lowest eigenvalues vanishes as h goes to zero. In particular, we have 
that: 



Lemma 3. Let 



A2 — Ai A2 + Ai 
LO = ana \l = 



ipR = ^ [ipi + ip2\ and ipL = ^ [ipi - (/32] 

where if 12 are the normalized eigenvectors associated to Ai 2. Then there exist two 
positive constants C and T , independent of h, such that 

WvrVlWoo < Cuj (7) 

and 

uj < Ce-^l^, Vfte (0,r]. (8) 

As a result if follows that 



and 



for some c > 0. 



lini w = (9) 

ri->o 



lim ^LJ™!i = c (10) 



Proof. In order to prove this Lemma we observe that is a symmetric double well 
potential with non-zero barrier between the wells. That is, let J > be small 
enough and let us define the two sets 

BR = {xeR+ : V{x) < Vrnrn + S} \ . ^r, ^ 

= |x e R : V{x) < Vmin +^} ) 
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From condition (H) it follows also that 

Br = [6, c] and Bl = [-c, -b] 
for some c > a > b > 0. The sets -B_r.l are usually called "wells". Let 

Ts = / ^/max[V{x) - (Kmn + 6),0]dx > 0, 

J-b 

be the Agmon distance between the two wells. From these facts and from standard 
WKB arguments (see § and §) then (@)— ® follow for some F e [To,Ts]. □ 

Remark 1. By definition it follows that ip]i{—x) = ipl{x); moreover, from (0). it 
follows that these functions are localized on only one of the "wells" Br and Bl, 
e.g.: 

/ \ipRix)\^dx = l+0{e-^/^) 

JBr 

for some C > 0. For such a reason we call them " single- well" (normalized) states. 

Remark 2. We underline that, by assuming some regularity properties on the 
potential V, then it is possible to obtain the precise asymptotic behavior of the 
sphtting as h goes to zero . 

2.3. Assumptions on the parameter. We assume that the two parameters uj 
and e are such that 

UJ ^ and e — > as S — > 
and there exists a positive constant C such that 

-<C, c^W^ll Vhe{0,h*]. (11) 

2.4. Assumption on the initial state. Let 

^c = I - {(PRr)'PR- {VL,-)V'L (12) 

be the projection operator onto the eigenspace orthogonal to the 2-dimensional 
eigenspace associated to the doublet {Ai_2}- Let -0° be the initial wavefunction, 
we assume that: 

Hypothesis 2. n^V'" = 0. 



3. Global existence of the solution and conservation laws 

Here, we prove that the Cauchy problem (H) admits a solution for all time 
provided that assumptions Hyp.l~2 are satisfied and the strength e of the nonlin- 
ear perturbation is small enough. Moreover, we prove a "priori" estimate of the 
solution -0. 

The following results hold. 

Theorem 1. There exist h* > and eo > such that for any h € (0,/i*] and 

e G [— eo,eo] then the Cauchy problem ^) admits a unique solution xp{t,x) G 
for any i G R. Moreover, the following conservation laws hold: 

Il^(i,-)I1 = I|0°(-)I1 = 1 (13) 
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5V' 



dx 



+ {V^J,^b) + -e\\i,Y 



(14) 



Proof. From Hyp. 2 it follows that 

From this fact and from Lemma 2 then ip'^ £ H^. Therefore, existence of the global 
solution ■0 € C(R, H^) and the conservation law s and (Q) follow from known 
results (see, e.g., the papers quoted in and ijlSj) for any e > (repulsive non- 
linear perturbation) and for any e G (— eo,0) for some eq > (attractive nonlinear 
perturbation) . □ 

Remark 3. There exists a positive constant C independent of h and e such that 

\E{4,)-Vr,,,n\<C{u + h + eh-^/^), Vne(0,r], Vee[-eo,eo]. (15) 

This estimate immediately follows from (p^, from the assumption Hyp. 2 and from 
Lemmas 1 and 2. Indeed, from Hyp. 2 it follows that 



0||4 
4 



where HV'^IU ^ Gh from (^) and where 

(77o(ci</7i + C2V52), (cit^Sl + C2V52)) = Ai | Ci ^ + A2 | C2 P 

From these facts and from (p^, then ( [T5| ) follows. 
Theorem 2. Lei eo(ft) &e a junction such that 

lim eo{h)/h^ = 0. 



r2 



2u;|c2| 



(16) 



The solution ip of equation (^ satisfies to the following uniform estimate: there 
exists a positive constant C independent of t, h and e, such that 



< C 



p-2 
4p 



Vp e [2,+oo], 



(17) 



dip 


< c 


dx 





\EilP°)-Vrmn\ 



h? 



for all time and Vft € (0, h*], Ve G [— eo(fi.), eo(?i)]. 
Proof. In order to prove the estimate ([l^) let 



/J = 



A 



i?(^°) - 



V2m' A:2 
Then, the conservation laws (O) and dl^) imply that 



dip 
dx 



1 



sign(e)]p2||V' 



2||2 



< A, 



where 



p=|e|l/2/^«l, 
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since (|l6|). In particular, if we set 

X = pip 

then the above equation takes the form 



dx 



dx 



i[sign(e)]||x'f <Ap2 



from which it foUows that 

2 



dx 



dx 



<p'm + \\\xY = p'm + \\\x\\\ 



(18) 



From the Gaghardo-Niremberg inequahty (see, for instance, [g[ and where the 
dimension is here equal to 1) 



|2(T + 2 ^ Q 



dx 



dx 



IIXll2<T+2 

where we choose cr = 1, it follows that 

dx 



Mr", Va>0, 



(19) 



11x111 <C 



dx 



\x\?<c 



dx 



dx 



P 



since 
that 



= /9 1 1 '01 1 = p and II?/' 1 1 = 1- By inserting this inequality in ( |18D it follows 
satisfies to 



dx 



dx 



<p2|A| + Cp3 



dx 



dx 



(20) 



for any t G R. From (g^) immediately follows that 
dx 



dx 



< V|A|p(l + o(l)) , as p->0. 



Hence, 



dx 



< 



C-\/|A| and, from (p^), we have that 



\mp<c 



9V 



dx 



a/p 



< C|A|(P^2)/*P 



where we choose now a — 2-^, i.e. p = 2a + 2. 



□ 



Remark 4. Condition ( p6| ) is true in the semiclassical limit and under the assump- 
tion (lul). 



Remark 5. From the fact £'(V'o) — Knm = 0{h), which follows from (||), ( [T5| ) and 
(H), and from the bounds (|l^) and (|ll|) then it follows that G R, V/i G (0, h*], 
Ve e [-eo(?i),eo(?i)] then 



|-0||p <Ch , ype [2, +oo], and 



dip 
dx 



(21) 
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4. Two-level approximation 

For our purposes it is more convenient to make the substitution ip e~*^*/''-0, 
hence equation (j^) takes the following form (where, with abuse of notation, we still 
denote the new function by ^p) 

ihiP = {Ho-n)i^ + e\ip\'^i^, ipix,0) = ^p°{x). (22) 

Let us write the solution of this equation in the form 

ip{t,x) = aR{t)ipii{x) +aL{t)ipL{x) +^c{t,x), (23) 

where ault) and ai(t) are unknown complex- valued functions depending on time 



and ipc — IIcV', He, defined in (|12|), is the projection onto the space orthogonal to 
the two-dimensional space spanned by the two "single- well" states cpn and ipL, i.e.: 

From the conservation law (|l^) it follows that 

\aR{t)\^ + |aL(t)p + WMt, -W = l,Vt e R. (24) 



By substituting -0 by ( |23| ) in equation (g) we obtain that a^j, ql and tpc must 
satisfy to the system of differential equations 

ihciL = -ujur + e{ipL, l^l'^ip) (25) 

zfiV-c = (Ho - n)^c + en.iv^pv- 

By substituting again tp by (^) in the first two equations of the above system 
then we obtain that these equations take the form 

ihcLR = -ujaL + ec\aR\'^aR + evR , , 

ihciL = -LoaR + ec\aL\^aL + erL 



where 

c^\\^lf = \\^ir = Oih-') (27) 
and where vr and are given by 

rR = {fR, l-iAlV) - \aR.\aR{fR, \(Pr\'^(Pr} 

= {fR, IV'IVl) + aj?(|</5_Rp, |0Lp + aRipR(j)L + aR(fR(j)L) 
rL = (y^L, I^AlV) - |ai|aL('/fL, |¥'l|V-l) 

= {ipL, \iP\'^4'r) + cilH^lI'^, + aL^L<t)R + aL<PL4>R) 



where 



6l = QLifiL + i'c and (pR = Oi^i^ii; -f V'c 



We denote two- level approximation the solutions bR and 6l of the system of 
ordinary differential equations 

ihhR = -ujbi + ec\bR\'^bR 

obtained by neglecting the remainder terms r/j and in (p6[). It is easy to see 
that the solution of this system satisfies to the following conservation law 

|6i?(i)l' + = \bRm^ + \hm? = |ai?(0)P + |aL(0)p = 1, (29) 
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and, moreover, it is also possible to explicitely compute (see and Appendix B 
in the solution of ( |2^ ) by means of elliptic functions. In particular, we obtain 
that the imbalance function, defined as 



is given by 



zit) = 



z{t)^\bn{t)\'-\bUt)f, 

Acn[Arj{ujt/h~ To)/2k,k], if fc < 1, 
Adn[Ari{ujt/h-To)/2,l/k], if k > 1, 



(30) 



where ry = ec/w, tq depends on the initial condition, 

' 77z2(0)/4, 



/ ^ Vl - z2(o)cos 
9 = arg(&ij) — arg(6i) is the relative phase, 

2V2 



A = 



1 , 
4^- 



1 + /r; - ( 1 + -/r; 



1/2 



and 



1 - 



l + i/7? 



(31) 



We underline that z{t) periodically assumes positive and negative values if, and 
only if, fc < 1. 

5. Stability of the two-level approximation 

Our main result consists in proving the stability of the two-level approximation 
when we restore the remainder terms and ri^ in equation (p8[). 
We prove that: 

Theorem 3. Let ipc ~ IIcV'; o,R{t) — {^i^r) and aL{t) — {t}),Lpi,), where -0 is the 
solution of equation ^^), let b^it) and b^it) be the solution of the system of ordi- 
nary differential equations (28). Let e G [— eo(/i), eo(fi)], where eo(?i) satisfies the 
condition (jJ^. Then, for any r' > there exists a positive constant C independent 
of e, h and t such that: 



\bR.,L{t) - aB.,L{t)\ <Ce 



-Ch" 



and Uc{-M < Ce 



-Ch- 



(32) 



for any h e (0, h*] and for any t G [0, hr' /uj] 



Proof. For the sake of simplicity, hereafter, let us drop out the parameters where 
this does not cause misunderstanding. In order to prove the theorem we introduce 
the " slow time" r = Lot/h and let 



Ar,l{t) = an^L{t) 

Br.,l{^) = bR^L{t) 

Then (|§) and 
respect to r) 



e ec 
Rr,l{t) = —rR^L[t) and rj = —. 

UJ UJ 



respectively take the form (here ' denotes the derivative with 

(33) 



A'r = 
A' = 



iAL - iTj\ARyAR + Rr 
iAR~ir]\AL\^AL+RL 
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and 

B'j^ = iBL~iv\BR\^BR 



(34) 



B'l = iBR-ir^lBL^BL 

satisfying to the same initial condition 

Br,l{0) = Ar^l{0) - afl,L(0), 

since ( p^ and (|2^), they are such that 

\BR{T)f + \BL{T)\' = l, \AR{T)f + \AL{r)f<l. (35) 

In a more concise way, with an obvious meaning of the notation, we can write ( |33| ) 
and (H) as 

A' = f{A) + R and B' = f{B), A{0) = 5(0) = a(0), (36) 
where A^BeS^ since (|3|), 5^ = {(zi,Z2) G : l^ip + |z2p < !}• 
Lemma 4. T/ie function f : satisfies to the Lipschitz condition: 

\f{A)^f{B)\<L\A-B\,L = l + 3Tj. (37) 
Proof. According with the notation (^ we have that 

|/(A)-/(i?)|=max[|/fi|,|/i|] 
where |A| < 1 and |-B| < 1 since A, B £ and where 

fR = {AL-BL)-vi\AR\^AR-\BR,\^BR) 
fL = {AR-BR)^r^{\AL\^AL-\BL\^BL) 



Then (37) immediately follows since 

fR = {Al ~Bl)~v [\Br\^{Ar ~ Br) + Ar{\Ar,\ + \Br\){\Ar\ - \Br\)] 

where \\Ar\ — \Br\\ < \Ar — Br\ and where the other term will be treated in 
the same way. □ 

Lemma 5. Let 

(3 — max[ce, lo] 



where c is defined in i\2l[). Let ipc = ^cip where ip is the solution of equation 
it satisfies to the following uniform estimate 

Uc\\ < C(3h-^'^ \eTi^[Ch-^/^{et/h)] + ij , € R, (38) 

for some positive constant C independent of h, e and t. 

Proof. As a first step we consider the following raw estimates: 

Uc\\p<Ch-^, Vpe[2,+oo], VieR, (39) 

and 

VrM <Ch-^'^, V<e R. 

Indeed, (^9|) immediately follows from the Minkowski inequality and from (pl|): 

cft-^ > ii^iip > - {\aRmwR\\p + WLmwdp) + iiv'cIIp 
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where |afl,L(i)l ^ 1 ^-i^^ where (^_r,l satisfy the bound (^. In the same way, from 
Lemma 2 and Theorem 2, it follows that 

\rR\ < C||^flV^2|| .||^||^||^^||4 

< C||^fl||oo||^||4llV^II+C||(^fl||| 

and similarly for \ri^\. 

Now, in order to prove the estimate ( ^8| ) we make use of the third equation of 
25|) from which it follows that 

since "ipc — ncV"" — from assumption Hyp. 2. 
Let ip = ip + Ipc where ip = aRipR + ul^l, then 



Therefore, we can write 



where 



/ 
// 
/// 



4'c = + II + HI] 



For what concerns the first term we have that, by integrating by part, 



"'0 



ds 



ds 



Let us underline that from Lemma 1 it follows that the following operators, from 
into L^, are bounded 



^-i{Ho-n){t-s)/h 



= 1, \\n[Ho-n]-^n,\\<c, 



and, from Lemma 2 and equations (|24| ) and (|2g) we have the following uniform 
estimate for any i e R 

_ 1 _ i _ 

ll^jllp < (|ar| + Iflil) (||'Pi?J|p + ||<Pl|Ip) < C/i" max[ce,u;,eh^'!]h 
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Then we have that 



|/|| < Cmax {||^3(^^.)||^^||^(^,^.)^2(^^.)||| 

se[Q,t] 

< C max {Ms,-)\\l + ms,-)\\-Ms,-)\\l} 
se[o,t\ 



For what concerns the other two terms we have that 

l|//||</ ||^c||-||mlloods<cft-i/2 / Wij^Wds 

JQ Jo 

since ||v//|loo < Ch^^^^, and similarly 

< / ll^cll • WfiiiWoods < Ch-^/^ f U,\\ds. 

Jo Jo 
Indeed, from Lemma 2 and (39) it follows that 

||'^//||oo < C{Ml + ||V^c||oo||<^||oo + ll^cIlL} < Ch~'^' 

and 

\\^iii\\oo<Ml<ch-^/'. 

Collecting all these results and denoting 

9{t)^\\M-M 

we have that g(t) is a positive real valued function satisfying the estimate 

, e 



g{t) < Cj\h~^'^ g{s)ds + h-^/^{l + th-^l3) 



< a g{s)ds + a + abt, a = C—j-, b^—. 
Jo n-^'^ n 

From this estimate, since V'c(O) = and from the Gronwall's Lemma (see 1^, pag. 
19) it follows that 

g{t) < a + abt + a [ e^^*"^) (a + abs)ds = -b + ae"* + be"* 
Jo 

< C/3 



^cem-^'^ + 1 



f,3/2 

proving so the result. □ 



From the inequality and from the assumption (|11|) it follows that for any 
fixed r' > then there exists C > satisfying the second inequality in (^). 

Lemma 6. For any fixed t' > then the remainder terms r^ and r^ satisfy to the 
following uniform estimate: 



max[\rR\,\rL\] < Cl3h-^e^''~"^\ Vt E [0,r'ft/ 
for some positive constant C independent of h, e and t. 



UJ\ 
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Proof. Let US only consider the term |ri^|, the other term \rL\ could be treated in 
the same way. By definition and since max[|ai^|, \aL\] < 1 then it follows that 

kfll < +|(^fl^L,|^P)| (40) 

+ |(^fl|VI',^c)| (41) 

+ Kl'/'Kp, |0Lp + aRipR^L + aR<pR(l>L)\ (42) 

and we estimate separately each term. 

From Lemma 3, equation ([l^ ) and the Holder inequality it follows that the term 
10|) satisfies to the following estimate: 

\{fR^L, < Wv'B^'fLWoc ■ \\lp'^\\l < Cuj. 

From Lemma 5 and the Holder inequality it follows that the term ( ^l|) satisfies 
to the following estimates: 

\{M^\\^c)\ < \\v^B.\\oo ■ m\ ■ \m < C/?rtcr.-v= 

and that the term (ji^) satisfies to the estimate 

|(|<^i?J^, |0Lp + aRLpR<f>L + aRipR<j)L) \ < 

< C [\\(pR(pL\\oa + ll'/'flJIooll'/'clP + ||'/'_R</'l||oo||'/'c||] < Cu 

Collecting all these estimates we obtain the proof of the Lemma. □ 



The proof of the Theorem is almost done. Indeed, equations (|3q ) can be 
rewritten in the integral form: 



A{r)^AiO)+ r f[A{s)]ds+ f 
Jo Jo 



Rds 



and 



B{t) - 5(0) + / nB{s)]ds 
Jo 

from which it follows that for any r £ [0, t'] 

\A{t)-B{t)\ < r\f[A{s)]--f[B{s)]\ds+ r\R\ds 



< a[ \A{s)~B{s)\ds + bT, a^L, b^C 
Jo 



since Lemmas 4 and 5. From this inequality and by means of the Gronwall's 
Lemma we finally obtain that 

\A(t)-B{t)\ < br + ab e"(*-'')sds = - [e°^ - 1] 

Jo a 

^ C eph-^e^^~"'" 
~ L u> 

proving so ( ^^ since 

" + '<i = l + 3ry<C'^, 



C't 

for some C" > 0, which implies that -/^ < C for some C > 0. □ 
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Remark 6. Recalling that uj = 0{e then the above theorem implies that for 
any a < I and for any r' > there exists C such that 

IHlW -afl,LWI < Cuj'^ and \\M;t)\\ < Cc^", Vt e [O^Ht'/uj]. 
6. Destruction of the beating motion for large nonlinearity 

6.1. The unperturbed case e = 0. Under the assumption Hyp. 2 it follows that 
the solution of the unperturbed equation 

is simply given by 



ip{t,x) = e 



cos[ujt/n) + I r= — sm(ujt/rL) 



V2 



V2 



+e 



-int/h 



cos[Lut/h) — I r= — sin(cjr/n,j 



V2 



where 



Cl,2 = (^l,2,^°), |C1P + |C2| = 



V2 



1. 



)t/h 



a periodic function with period 



Hence, ■(/'(t, x) is, up to the phase factor e~**^^" 
T = T:h/uj. 

In particular, if -0 initially coincides with a single-well state, e.g. — Vb.^ then 



and the state x) performs a beating motion. That is the state, initially localized 
on the well Br, is localized on the other well Bl after half a period and, after a 
whole period, it "returns" on the initial well, and so on. In particular, let us 
consider the motion of the "center of mass" defined here as 



X(x)\^l:{t,x)\^dx 



where X G C[R) n -/^^(R) is a given bounded function such that X{—x) ~ —X{x). 
We have that 

{Xy ^[cos^{ojt/h)-fim^{ojt/h)\ I X{x)\ipB,{x)\'^dx 

Jr 

is a periodic function which periodically assumes positive and negative values, i.e. 
we have the well know beating motion for double-well problem. 

6.2. The perturbed case e ^ 0. In such a case it follows that the "center of 
mass" is given by 

{Xy = c[|afl(i)|' - |aL(t)P] +r, c= (XipR, ^r), 

where the remainder term r satisfies the uniform estimate 

|r| = 2\lk[aRdL{X^R,ipL) + {X^,i^c)]\ 

< 2[||^fl^i||,o + ||^||oo||V'IIIIV'c||] 

< Ce-^'''\ yt e [Q,nT'/uj]. 

If we denote by z{t) the imbalance function defined in (50) then it follows that 



\aRit)\^ -\aL{t)\^ ^z{t) + Ce 



-Ch- 
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hence 

{Xy ^ cz{t) + Ce~^^'\ yte [0,hT'/uj]. 
Then, we have that: 

Theorem 4. Let Hyp. 1 and 2 he satisfied. Let he defined as in (pip^ it depends 
on the initial wavef unction ip^ . Let r' > fixed, (X)* is, up to a remainder term, 
a periodic function for any t G [0, St'/w]. In particular, if: 

i) fc^ < 1 then (X)* periodically assumes positive and negative values (i.e. the 
heating motion .still persists); 

ii) > 1 then (Xy has a definite sign (i.e. the beating motion is forbidden). 

Remark 7. Let us close by underhning that when the wavefunction is initially 
prepared on just one well, e.g. = 'PRi then 

I^~\ri and k^ ^ ^if . 

Therefore, from the theorem above it follows that for I77I larger than the critical 
value 4 the beating motion is forbidden. In such a way, we put on a full rigorous 
basis the results obtained by in the two-level approximation. 
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